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Polarization in the inner region of Pulsar Wind Nebulae 
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Abstract. We present here the first effort to compute synthetic synchrotron polarization maps of Pulsar Wind Nebulae (PWNe). 
Our goal is to highlight how polarization can be used as an additional diagnostic tool for the flow structure in the inner regions 
of these nebulae. Recent numerical simulations suggest the presence of flow velocities ~ 0.5 c in the surroundings of the 
termination shock, where most of the high energy emission comes from. We construct polarization maps taking into account 
relativistic effects like Doppler boosting and position angle swing. The effect of different bulk velocities is clarified with the 
help of a toy-model consisting of a uniformly emitting torus. We also present a map based on recent numerical simulations 
of the entire nebula and compare it with presently available data. The comparison with upcoming high resolution observations 
could provide new insight into the inner structure of the nebula and put constraints on the geometrical properties of the magnetic 
field. 
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1. Introduction 

Pulsar Wind Nebulae (PWNe) are bubbles of relativistic par- 
ticles and magnetic field produced by the interaction of the 
relativistic pulsar wind with the ambient medium (ISM or 
Supernova Remnant). These bubbles are souces of non thermal 
(synchrotron and Inverse Compton) emission, in a broad range 
of frequencies from radio wavelengths to gamma rays. X-ray 
images of PWNe show an axisymmetric feature known as the 
jet-torus structure. This feature have been observed by now in 
a number of PWNe, among which the Crab Nebula (Weisskopf 
et al. I2000L Ve la (He lfand et al. I2QQTI and PSR B 1509-58 
(Gaensler et al. 120021) . Recent relativistic MHD simulations 
(Komissarov & Lyubarsky 2003, Komissarov & Lyubarsky 
120041 Del Zanna et al. 2004) have shown that this structure 
naturally arises if the pulsar wind energy flux is anisotropic, 
higher in the equatorial plane of the pulsar rotation than along 
the polar axis, and the ratio of Poynting to particle energy flux, 
<x, is greater than approximately 1%. We will not discuss here 
the details of those simulations and the related flow structure 
inside PWNe. Let us just recall that the derived synthetic syn- 
chrotron maps allow to recover, at least qualitatively, the torus 
structure with a bright main arc and other features, like the knot 
that is observed in the Crab Nebula almost coincident with the 
pulsar position. 

Despite the success of these models at reproducing most of 
the basic observational features, there are still details that can- 



not be accounted for because of the intrinsic model assump- 
tions. A major limitation appears to be the hypothesis of a 
purely toroidal field, which vanishes on the axis, thus strongly 
suppressing the jet emissivity. In order to reproduce the ob- 
served PWN jets an additional magnetic field component, ei- 
ther poloidal or disordered, must be necessarily introduced. At 
present, in the absence of fully 2.5-D, or even 3-D, MHD sim- 
ulations, only ad-hoc prescriptions have been used. Here we 
look for constraints that can be put on the geometry of the mag- 
netic field, in particular the presence and strength of a poloidal 
component. 

As shown by Shibata et al. 120031 a completely disordered 
(isotropic) magnetic field has a very different emission signa- 
ture with respect to a purely toroidal one, in the case of low bulk 
speed in the torus. However, the flow speed in the torus region, 
both from observations and the above cited MHD simulations, 
is of the order of ~ 0.5 c and Doppler boosting effects strongly 
affect the overall emission, making it difficult to constrain the 
field structure. 

On the other hand, polarization is supposed to provide a 
more direct way to investigate the magnetic field configura- 
tion. Radio polarization maps are available for many PWNe, 
among which the Crab Nebula (Wilson lT972l Velusamv l 19851) . 
3C58 (Wilson & Weiler [T9761 W eiler & Seiesland IT97TI , 
G328.4+0.2 (Gaensler et al. l2000l) . Vela ( Dodson et al. l2003l . 
SNR G320.4-01.2 (Gaensler et al. but radio emitting 

particles have too long synchrotron cooling times, so that the 
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nebular emission is usually dominated by the outer regions (the 
Crab Nebula torus is barely visible in radio). Radio polarization 
might then just provide a good estimate of the degree of or- 
dered versus disordered magnetic field for the overall nebula, 
but cannot be used to constrain the region close to the termi- 
nation shock. What one would need is a polarization map at 
higher frequencies. In the X-rays, where the jet-torus structure 
was first clearly recognized, polarimetry is unfortunately not 
available at present. However, the torus and other bright fea- 
tures near the termination shock (like the inner ring, the time- 
dependent wisps, and the central knot) are well visible in the 
Crab Nebula at optical wavelengths, at which high resolution 
polarimetry might be employed. Until now only few optical 
polarization maps of the Crab Nebula (Schmidt et al. 119791 



Hickson & van der Bergh 1990, Michel et al. 1991 ) have have 
been published, in general with too low spatial resolution to in- 
vestigate the inner region in detail, though recent Hubble data 
appear to be more promising (James Graham, personal com- 
munication). 

Here we present synthetic polarization maps as a novel di- 
agnostic tool for the physics of PWNe. We will first discuss the 
basic ingredients that determine the polarization direction and 
fraction within the framework of a simplified toy-model. Then 
this tool will be applied to the flow properties derived from (ax- 
isymmetric) relativistic MHD simulations, as those presented 
in Del Zanna et al. 2004| Our work is aimed at clarifying if 
and how polarization maps might be used to infer information 
on the inner flow structure, by comparison with high-resolution 
optical observations. 

In Sect. 2 we describe the model adopted for the emission 
and for the polarization, both for the toy-model and in the case 
of the more complex flow pattern derived from numerical sim- 
ulations. The results are discussed and the effect of relativis- 
tic bulk flow speeds on the polarization are shown. Finally in 
Sect. 3 we summarize our conclusions. 

2. Synchrotron emissivity and polarization 

Let us briefly recall here the derivation of the local synchrotron 
emissivity s(v), at a given frequency v, for the case in which the 
particles' energy distribution function at the termination shock 
is a power law 



f(e ) oc P 6 



,-(2or+l) 



(1) 



where Po is the total thermal pressure immediately downstream 
of the shock. As shown in Bucciantini et al. 2005 (to which 
the reader is referred for a detailed derivation of the following 
equations), for a close to 0.5 the emissivity at a generic position 
inside the nebula becomes 



s(v) oc PD a+z (B'J 



\a+ly-a 



(2) 



where P is the local thermal pressure, B' ± the magnetic field 
component perpendicular to the line of sight as measured in 
the fluid frame, and D is the Doppler boosting factor. 

Since we are interested to high energy emitting particles, 
whose characteristic lifetime is comparable to or shorter than 
the age of the nebula, synchrotron losses must be taken into 



account. Hence we correct the emissivity by multiplying the 
Eq.|2|by: 



[1 - €(y)/e 00 ] 
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where e(v) is the energy of particles emitting at frequency v 
(averaged over the pitch angle for an isotropic distribution) 
and 6oo is the local energy of a particle injected at the termi- 
nation shock with infinite energy (we assume zero emissivity 
where 6oo < e(v)). The nebular emission map is then finally 
constructed by integrating along the line of sight. 

In order to build a polarization map, given a pattern of ve- 
locity streamlines, one also needs to evaluate the Stokes param- 
eters U, 2, V. However, since synchrotron radiation from high 
Lorentz factor particles is known to have a very low degree of 
circular polarization, here we will neglect it completely (V = 0) 
and thus evaluate only U and Q. 

By naming dl the total intensity of the radiation emitted by 
a fluid element (the emissivity s(v) in Eq. 0), the value of dQ 
and dU are defined in terms of dl and the polarization angles 
as: 



dQ- 



a + 1 
a + 5/3 



dl cos 2%, 



dU= " 4 d/sin2r. 
a + 5/3 



(4) 
(5) 



We define x as the angle, in the plane of the sky, between the 
emitted electric field and the projected axis of the nebula. Its 
value is computed in two steps, as in Pariev et al. 2003] (see 
also B landford & Konigl [T979l Biornsson [T982l Qian & Zhang 
2003 ). First we evaluate the inclination angle if/ of the emitted 
radiation electric field e with respect to a plane containing the 
observer direction n and the velocity vector ft (considering that 
in the comoving frame e is normal both to the local magnetic 
field and to n). By choosing a Cartesian coordinate system with 
the x axis in the direction of fi and y parallel to n x ft we find 



tan iJj = 



B z ]fi\-fi-n/\fi\ 
B y 1-jffn ' 



(6) 



where we have assumed /?B = 0=># X = 0, since our simpli- 
fied model and the MHD simulations both assume a geometry 
with a purely toroidal field and a purely poloidal velocity. Note 
that the above equation differs from Eq. (24) in Qian & Zhang 
120031 because of a mistake in the definition of the magnetic 
field components. Finally^ is obtained by adding to if/ a quan- 
tity representing the angle between the plane containing ft and 
n, and the plane containing n and the projected axis of the neb- 
ula. Once the value of x is known for the radiation emitted by 
a fluid element, the Stokes parameters are computed by inte- 
grating the equations for dQ and dU along the line of sight, 
thus providing maps of both polarization fraction and position 
angle. 

2.1. Uniform torus model 

For the sake of clarity let us first consider the simple case 
where the emission is concentrated in a homogeneously radiat- 
ing torus around the pulsar. The existence of arcs and ring-like 
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features observed in a number of PWNe may be modeled as a 
first approximation in this way. In a spherical reference system 
(r, 0, (p) centered on the pulsar position the torus can be shaped 
as: 

(r sin - R x f + (r cos 6) 2 -R 2 2 <0, (7) 

where R\ is the torus principal radius and R2 is the radius of the 
cross section (both will be given here in arbitrary units). Inside 
this region we consider quantities with homogeneous values, 
with the magnetic field assumed to be purely toroidal, whereas 
the flow velocity is assumed purely radial. In this simplified 
model we do not take into account the effect of synchrotron 
losses, that is the correction in Eq. 0, thus the resulting maps 
are for a generic frequency. Outside the torus we impose zero 
emissivity. 

Concerning the polarization fraction we find that in the en- 
tire torus it is very close to its theoretical maximum. When the 
bulk flow exceeds 0.2 c, Doppler boosting effects are impor- 
tant and even where contributions from the front and back side 
of the torus are integrated together along the line of sight the 
depolarization is still negligible there. Minor deviations form 
the maximum polarization fraction are visible just on the outer 
sides. 

More interesting is the behavior of the polarization vector. 
As expected, close to the symmetry axis of the torus, coinci- 
dent with that of the nebula, this is aligned with the axis itself 
(due to the assumption of a purely toroidal field), and the incli- 
nation of the polarization vector with respect to the symmetry 
axis increases moving away from it. The position where the 
polarization vector is perpendicular to the axis is a function of 
both the inclination of the torus with respect to the plane of the 
sky 6, and the flow speed. 

We consider specifically an inclination of 6 - 30°, cor- 
responding approximately to the best estimate for the Crab 
Nebula. We find that for increasing values of the bulk speed 
of 0.2, 0.4, 0.6, 0.8 c, the corresponding position at which the 
polarization vector is perpendicular to the axis moves closer in 
(see fig.[D- We estimate the position of this transition in a polar 
coordinate system centered on the pulsar and lying on the plane 
of the sky, with angles measured from the axis of the nebula. 
For the above mentioned values of velocity, the angle of the 
transition is 90°, 80°, 70°, 60°, respectively. 

Looking at Fig. 5 of Schmidt et al. ( 1979) we can notice 
in the regions labeled "Wl" and "W2" that the electric field is 
perpendicular to the nebular axis at an angle of about ~ 70°. 

Besides, the polarization vector in the back side of the torus 
remains aligned with the axis, while diverges from it in the 
front side. This effect increases with bulk flow speed. Again, 
such difference in the divergence of the polarization vector be- 
tween the front and back side of the nebula is observed in Fig. 5 
of Schmidt et al. 119791 If one wants to reproduce the observed 
front-back difference in the divergence of the electric field vec- 
tor, a bulk velocity greater than at least 0.3 c is required. We 
thus tentatively conclude that evidence for high flow velocities 
in the nebula was already present in that map. 



2.2. Numerical simulations 

Synthetic polarization maps to compare with the data can be 
obtained by applying the scheme outlined above to more so- 
phisticated models for the magnetic field structure, flow ve- 
locity pattern, and emitting particle distribution. Axisymmetric 
relativistic MHD simulations (Del Zanna et al. 2004) provide 
the full set of dynamical quantities. In order to properly eval- 
uate the emission at high frequencies, as discussed in the pre- 
vious section, what is left to take into account is the effect of 
synchrotron losses. In order to do this we have added to our 
code a tracer that evolves €00 (see Eq. Q) from the termination 
shock to the outer regions of the nebula following the stream- 
lines, thus allowing us to model the modifications in the emit- 
ting particles' distribution function. 

We show in Fig. [2] the results for the total polarization frac- 
tion and the polarization vector direction, at optical frequen- 
cies. Values are computed based on a recent simulation per- 
formed using a fixed spherical boundary condition (imposing a 
constant, non relativistic, flow speed Vpwn at the nebular outer 
boundary), thus neglecting the interaction of the PWN with the 
expanding ejecta. In the simulation we have used the following 
wind parameters: anisotropy a = 0.1, magnetization <x = 0.03, 
and b = 10 (b parameterizes the width of the striped wind 
region of reduced magnetic field around the equator, see Del 
Zanna et al. 2004 for a precise definition of these parameters). 
The other quantities are typical of the Crab Nebula: luminos- 
ity L = 5 x 10 38 erg s" 1 , age t = 1000 yr, boundary speed 
Vpwn = 1500 km s" 1 . 

Let us discuss the results in Fig. [2] and their implications. 
First of all we identify two main structures in the optical sur- 
face brightness map (right panel): an inner ring and an outer 
torus. As in the uniform ring case, the polarization reaches its 
maximum on the axis. The front side of the inner ring has a 
very high degree of polarization in a large region, up to about 
60° from the symmetry axis of the nebula, while, in the back 
side, the polarization tends to decline faster moving away from 
the axis. We also observe an enhancement of polarization in 
the region corresponding to the torus where values are of order 
70% of the maximum. Contrary to the toy-model case, where 
we found the entire torus to be highly polarized, now the polar- 
ization fraction drops below 50% along the sides of the inner 
ring, and there are regions at about 80° where the emission is 
almost depolarized. 

Concerning the direction of the polarization vector, we 
again observe the effect of Doppler boosting in the inner ring, 
where the polarization vector becomes perpendicular to the 
axis of the nebula at about 70°. Taking as reference our toy- 
model, this would imply a speed ~ 0.6 c. The same effect is 
still present but less strong in the torus, where the bulk flow 
speed is much lower. It also seems that in the back side of the 
ring the polarization tends to remain parallel to the axis up to 
larger angular distances than in the front side. 

In the outer region of the nebula the polarization vector is 
in the radial direction, as expected for a slow moving flow em- 
bedded in a toroidal field. However, the results in this region 
should not be taken too seriously because here our neglect of 
the interaction with the SNR material may be crucial. 
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Fig. 1. Polarization vector and total emission (normalized to the maximum) in the uniform torus case (here R\ = 12.5 and R2 - 5, 
in arbitrary units, chosen to reproduce approximately the ratio ^2/^1 estimated for the Crab Nebula torus). The inclination angle 
of the symmetry axis with respect to the plane of the sky is S = 30°. Top-left v = 0.2 c\ top-right v = 0.4 c; bottom-left v = 0.6 c; 
bottom-right v = 0.8 c. The length of the polarization vector is proportional to the polarization fraction. 



Our results suggest that the behavior of the polarization 
vector in the inner ring (i.e. the swing in position angle) might 
be used to constrain the flow speed, while the polarization frac- 
tion may help estimating the amount of disordered or poloidal 
magnetic field. 



3. Conclusions 

In this paper we have presented preliminary results concerning 
the properties of non-thermal emission of PWNe in 2-D also 
computing, for the first time the (linear) polarization structure. 
The latter task, on which this paper is mainly focused, has been 
carried out both in the framework of a simplified model and by 
using the results of relativistic MHD simulations of the flow 
inside the PWN (Del Zanna et al. 2004). Our results suggest 



N. Bucciantini, L. Del Zanna, E. Amato, D. Volpi: Polarization in the inner region of Pulsar Wind Nebulae 



5 



Polarization Fraction 



Polarization Diraction 





R ■ WBKm—. . . . . -dd 




-1.5 












-1.5 


-1 .0 


-0.5 0.0 0.5 
x (ly) 


1 .0 1 


5 


-1.5 


-1.0 


-0.5 


0.0 

x (ly) 


0.5 


1 .0 1 


















0.0 


0.2 


0.4 0.6 


0.8 1 







0.2 


0.4 


0.6 




0.8 1 



Fig. 2. Results from full relativistic MHD simulations. Left panel: polarization fraction in the inner region of a PWN, normalized 
to the theoretical maximum (~ 70% for a particle distribution index a = 0.5). A: main torus; B: inner ring. Right panel: 
total optical emission (normalized to the maximum) and polarization vectors. The length of the vector is proportional to the 
polarization fraction. This map accounts for synchrotron losses, using typical values of age and size for the Crab Nebula. 



that the bulk flow in the emitting region and the polarization 
properties are strongly correlated. 

Our simplified model consists of a homogeneous torus with 
prescribed shape and inclination with respect to the plane of 
the sky. The torus has both uniform flow velocity, assumed 
to be purely radial, and magnetic field, purely toroidal. Under 
these conditions, we find that the polarization fraction is always 
close to the expected maximum, regardless of the radial speed, 
whereas, concerning the position angle, we have shown that for 
high velocities, due to relativistic effects, the polarization vec- 
tors diverge rapidly from the axis of the nebula, in the front 
side, but tend to remain aligned in the back side. 

A more complex pattern is found in the polarization 
maps based on numerical simulations. The synchrotron sur- 
face brightness map at optical frequencies shows two ring-like 
structures: a bright inner ring, at about the distance of the wind 
termination shock, and a fainter torus farther away from the 
pulsar. First of all we find a correlation between the intensity 
and the degree of polarization. This is however generally lower 
than in the previous case, with the emission almost completely 
depolarized at the edges of these brighter structures. As far as 
position angle is concerned, polarization swing is stronger in 
the inner ring than in the torus, due to a higher flow speed in 
the former, while both retain the front/back asymmetry found 
in the toy-model. 

We would like to point out that this asymmetry was already 
showing in the observations of the Crab Nebula, and could have 
been used to infer relativistic speeds in the inner parts of the 
nebula. It should be clear at this point that the combination of 
optical (and hopefully X-ray) high-resolution polarimetry, and 



synthetic polarization maps from relativistic MHD simulations, 
is promising to provide a powerful diagnostic tool for the flow 
structure in PWNe. Polarization measures might also help to 
clarify the origin of some puzzling emission features like the 
knot of the Crab Nebula, revealing whether they are associated 
with the nebular post-shock flow or a result of some kind of 
plasma instability in the pulsar wind. 

In order to understand how emission features and polariza- 
tion properties can be used to infer the wind structure and mag- 
netization, a more detailed study based on a set of numerical 
simulations is currently under progress. Results will be pre- 
sented in a future paper. 
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